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ABSTRACT 



Anelastic light scattering is computed numerically for model disordered sys- 
tems (linear chains and 2-dimensional site and bond percolators), with and 
without electrical disorder. A detailed analysis of the vibrational modes and of 
their Raman activity evidences that two extreme mechanisms for scattering may 
be singled out. One of these resembles scattering from finite size systems, while 
the other mechanisms originates from spatial fluctuations of the polarizability 
and is such that modes in even small frequency intervals may have very different 
Raman activities. As a consequence, the average coupling coefficient C{uj) is 
the variance of a zero-average quantity. Our analysis shows that for both linear 
chains and 2-dimensional percolators the second mechanism dominates over the 
first, and therefore Raman scattering from disordered systems is essentially due 
to spatial fluctuations. 
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1. INTRODUCTION 



Inelastic light scattering is a powerful tool to obtain information on the 
vibrational properties of solids. In the case of crystals, the existence of transla- 
tional symmetry produces a selection rule which links the pseudo-momentum of 
the created or annihilated phonon, k, to the momentum q exchanged between 
the photon and the system, i.e. k — q; a,s a consequence acoustic phonons pro- 
duce Brillouin scattering, while optical phonons at the centre of the Brillouin 
zone (fc <C 7r/a where a is the interatomic distance) produce Raman scattering. 
The reason why phonons with k ^q do not scatter, is that in this case the elec- 
tromagnetic waves produced by the polarizability modulations that are induced 
at different sites of the crystal by the modes, interfere in an exactly destructive 
way. 

In disordered solids the situation is quite different. Brillouin scattering still 
produces a sharp and intense peak as in crystals. In fact, the wavelength of the 
light A ~ 5OOOA is much greater than the correlation length of disorder so that 
the vibrational modes with such wavelength are much the same as phonons in 
crystals, i.e. plane waves which propagate in the disordered medium with a well 
defined fc. So the rule fc = g is still valid. In the case of Brillouin scattering 
it is customary to speak about "coherent" scattering because the polarizability 
changes induced by the vibrations with k ^q scatter the light inelastically in a 
coherent way, which means that the electric field amplitudes scattered by atoms 
placed at different sites sum up in phase. In the following we will assume g ~ 
and we will not consider the Brillouin scattering any more in this paper. 

Besides the Brillouin peaks, a continuum spectrum is present; this is due 
to the absence of spatial periodicity. In fact, the spatial fluctuations of the 
acousto optical constants prevent the above mentioned completely destructive 
interference among the scattered waves; this produces inelastic scattering which 
originates from the spatial variations of the electric polarizability (disorder- 
induced light scattering, DILS) and which has a continuum spectrum ranging 
from cj = up to some hundred wavcnumbcrs. This is often called "incoherent" 
scattering, though in our opinion this terminology can produce confusion. 

In the harmonic approximation the DILS scattered intensity I{uj,T) can be 
put in the form 

I{uj) cx [n{ijj) + 1]C{lo)p{uj)/uj 

where n{Lo,T) is the Bose-Einstein population factor, p{uj) is the density of 
vibrational states and C{uj) is the average light-vibration coupling coefficient 
for the modes whose frequency lies between oj and lo -\- Slu: 



where 



EE 



Q/3 
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BM\' (2) 



which defines the quantity Baf}{p) whose utihty will become clear later. Note 
that since Bap{p) is a linear function of the atomic displacements e{j\p)/ ^JrfTj 
produced by mode and since the sign of the e(j|p)'s is arbitrary, the only 
physically relevant quantity is \Bai3(p)\- 

For low frequency acoustic modes (i.e. those whose wavelength is greater 
than the correlation length of disorder) what is expected is C{uj) oc tj^ which 
has been observed experimentally in some systems Q . At higher frequency the 
situation is more complicated. 

The DILS Raman scattering is due to a spatially disordered distribution of 
the polarizabilities of the scattering units (ions, atoms or molecules). In infinite, 
mechanically and electrically ordered systems the scattering intensity vanishes. 
In disordered systems the spatial polarizability fluctuations produce a finite 
scattering intensity: these fluctuations can be due to mechanical disorder alone, 
to electrical disorder alone, or to both. The situation is somewhat different in 
finite systems, whose edges can be considered as regions where the polarizability 
varies from the value of the scattering units to that of vacuum, so that even fully 
ordered finite systems scatter. 

Even though scattering arises in all cases from polarizability variations, it 
may be useful to define two extreme types of scattering, (i) Scattering that 
is produced by spatial polarizability fluctuations, in the sense that different 
sites have different polarizabilities. This kind of scattering is present both in 
mechanically ordered Q] and disordered solids, {ii) "Finite-size" scattering, i.e. 
the only contribution to scattering present in fully ordered finite systems; this 
is due to edge effects. 

As we will see, even in mechanically disordered systems some of the scat- 
tering intensity can be ascribed to mechanisms very similar to edge effects. We 
will also see that these two scattering mechanisms have different characteristics 
which can be used to define them a posteriori. 

In the present paper we will study the relative importance of these two 
mechanisms; this depends in principle on the scattering system and on the 
frequency range. We will investigate two simple model systems: linear chains 
and two-dimensional site- and bond-percolators: the statistical method set up 
for linear chains will be utilized to study the percolators, which, at threshold, 
are fractal structures. 

Raman scattering from supposedly fractal systems was studied in a number 
of recent experimental |^ |^ papers. The experimental results seem to 
imply that the Raman coupling coefficient follows a power law C{oj) oc w^. 
Much numerical work [ p||pt [10| was devoted to determine x in percolators, while 
11 |12|| 13[ were devoted to derive scaling laws for C{uj) in 



theoretical efforts 
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terms of the mechanical parameters (hke the fractal dimension Z?, the spectral 
dimension d, and others) In the next sections, we will show that for 

the systems under consideration the DILS mechanism is by far the dominant 
one. 

The paper is organized as follows: in section 2 we present the results relative 
to linear chains with electrical and/or mechanical disorder. In section 3 we 
report the results on percolators, while the final section is devoted to discussion 
and conclusions. 



In order to better illustrate the scattering mechanisms proposed in the In- 
troduction, we will apply the previous formalism to various situations relative 
to the linear chain. 

Consider an ordered, infinite {N — > oo) linear chain with equal masses and 
equilibrium positions = aj, j — I....N, connected by equal springs. In this 
case, disregarding Brillouin scattering, and omitting the polarization subscripts, 
equation (g) becomes: 



where Aij = dir'^/du^ . The explicit form of Aij depends on the scattering 
mechanism used (Dipole-Induced-Dipole (DID), Bond Polarizability (BP), etc) 
but in any case it is an odd function, Aij = —Aji. 

It is easy to see that B{p) — Q due to cancellation effects which involve all 
the scatterers. In the case of a finite chain, on the contrary, B{p) ^ due 
to the absence of complete cancellation for the scatterers at the edges. It is 
clear that ordered systems do not scatter, and it is only the breakdown of the 
translational symmetry of the polarizability which makes finite ordered systems 
slightly Raman active. 

Another interesting case is that of electrical disorder in a mechanically or- 
dered system; the coupling coefficient was derived in here we study the 
distribution B{p) for one dimensional systems. Let us consider a linear chain 
of identical masses linked by identical next-neighbour springs of rest length a, 
and let us assume that the scattering mechanism is next-neighbour DID (ND) 
[|j (extension to three dimensions and full DID (FD) is straightforward). The 
electrical disorder is produced by randomly assigning to each mass a different 
isotropic, point-like polarizability a^. In this case 



2. LINEAR CHAIN 




2 




A, 




a'a'T^^\a) if i=j+l 
-a*aJr(3)(a) if i=j-l 



otherwise 
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where 

= -V„V/5V^(l/|r-|) (4) 

It is easy to obtain, with straightforward but lengthy algebra, the following 
expression: 



B(p) = ^T(3)(a)[e^'^''° - 1] 



(l + e-''=''°)(a)i^ + G 



(5) 



where 

F = ^foje*^^"^ (6) 
j 

G = ^(5aj(5aj+ie*'^''"^ (7) 
j 

For a large number of masses with polarizability and for many independent 
statistical replicas, the statistical distribution of the quantity is a gaussian 
centered at (F) = and with a variance cr|, = (F^) ~ N{{da)'^). The distribu- 
tion of G is also a zero-centered gaussian with CTq — {ap)'^/N. Therefore the 
distribution of B{p), which is the convolution of two gaussians, is again a zero 
centered gaussian with a variance cr^ given by: 



4 = 



N 

8 



= 4^(^j [(">'((H'>^*n'(fcpa) + {{5a)Ysin\k,al2)] (8) 

in agreement with eqs. (22), (31), (35) and (36) of reference [Q. 

From this analysis it is clear that the origin of DILS is in fluctuations: in 
fact, C{p) is just the variance of the zero-average quantity B(p). 

We consider next a linear chain with mi = 1 and m2 = 2 randomly dis- 
tributed with equal probability and linked by identical springs of elastic con- 
stant K = \, with periodic boundary conditions. The maximum frequency is 
^max = 2^Kl [I where [i is the lighter mass, and in the following lo will be given 
in units of LOmax- In the infinite chain all modes would be localized but in 
a finite chain the localization length of the lowest energy modes may be larger 
than the chain length, so that these modes appear not to be localized. Roughly 
speaking, there are two kinds of modes: 

• acoustic-like modes; 

• very localized modes, 

We assign the polarizabilities ai , a-i to the masses mi and respectively. 
Vibrational eigenvalues and eigenvectors were obtained by diagonalizing the 
dynamical matrix of many replicas of the system. In Fig. 1(a) we report Z3(|i?|), 
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i.e. the distribution of the |-B(p)| relative to the 2 modes with lowest frequency 
of 2000 chains containing 150 masses each, with ai = a2 = 2. C{uj) may be 
then obtained from equation (|^) by averaging on the modes with energy u!±Slo 
which are produced by the different realizations. 

In order to test if D{\B\) is a gaussian, we computed its first 10 moments: 



and compared them with the corresponding ones, M^, of the gaussian having 
the same zero and second moment. The ratios Mn/M^ are reported in Table 
I together with the test values. We find that the gaussian distribution is a 
very good approximation for all acoustic-like modes (i.e. for u < 0.5), so that 
the introduction of mechanical disorder does not change the results we obtained 
for electrical disorder alone. For some frequency intervals for tu > 0.5 important 
departures from gaussian shape are observed which show up in the appearance 
of peaks at well defined \B\ values (see Fig. 1(b)). 

By choosing narrower frequency intervals on which the averages are taken 
it is possible to enhance the peaks (see Fig. 1(c)). The frequencies around 
which Fig. 1(a) and Fig. l(b-c) were computed are indicated by the arrows in 
Fig. 2, where C(a;) is shown. The modes which produce the peaks at \B\ ^ 
are those where two neighboring light masses are surrounded by two long 
sequences of heavy masses, and oscillate in counter phase with ld k, 0.9, while the 
heavy masses are practically stationary. These modes scatter by a mechanism 
analogous to that of finite size systems: they cannot produce exact cancellation 
of the scattering amplitude. 

The present analysis evidences therefore that there are two extreme mech- 
anisms for Raman activity in disordered linear chains: DILS and finite-size 
scattering. The former produces a distribution D(\B\) which is a gaussian cen- 
tered at zero, while the latter produces peaks. In any case we have seen that the 
DILS mechanism is the largely dominant one when many masses are involved in 
the mode. Only in very narrow and high frequency ranges one can find modes 
localized on few atoms which give finite-size scattering. 

3. SCATTERING FROM 2-DIMENSIONAL PERCOLATING 

STRUCTURES 

The same analysis as in the previous section was performed on 2-dimensional 
site and bond percolators at percolation threshold, containing identical masses 
with identical bare polarizability. In the systems, each mass is bound to its 
nearest neighbours, the connecting springs are identical and scalar elasticity is 
assumed. Periodic boundary conditions are imposed and the eigenvalues are 
normalized to the maximum frequency of a full square lattice lOmax = 2-\/2i^/ \i. 
The Raman coupling coefficient was computed assuming that the scattering 




n 



6 



mechanism is full DID, nearest neighbour DID, and Bond Polarizability (BP) 

IE- 

In two dimensions equation ^ for FD becomes 



ij 7 



2 

2 



= \Bo.(sipW (9) 



To calculate the depolarized Raman scattering, B^yij)) (or Byx{p)) are used, 
while for the polarized scattering Bxx{p), B„j{p) and Byy{p) are used. We 
have computed separately the statistical distributions for the two groups of i3's 
(Bqq, or Baj3 with a ^ (3). In Fig. 3 (a-c) we report the B^a distributions 
relative to the modes in three frequency intervals (0.08-0.12, 0.37-0.43, 0.57- 

0. 63, respectively) for 450 rephcas of a 20x20 site percolator. For the two 
lowest-frequency intervals the distributions are very well fitted by gaussians, as 
shown in Table I. For the frequency range 0.57 < a; < 0.63, important deviations 
from the gaussian shape are observed near the origin, and are quantified in Table 

1. In particular, there is an excess in the first channel, showing that many modes 
are little Raman active, which means that cancellation effects are very relevant. 
By expanding the abscissa (Fig. 4) it comes out that most of the excess modes 
are confined to a very small range of values, of the order of w 10^^ times 
the value of a of the distribution of Fig. 3 (c) . 

In order to gain more insight on the nature of these low activity modes, 
we tried to characterize their wavefunction by the localization length evaluated 
according to the Thouless definition: 

where D is the fractal dimension and R is the participation ratio: 

In Fig. 5(a) we report the distribution of I for the « 10"* fractons with frequency 
0.57 < u! < 0.63, whose D{\B\) is shown in Fig. 3(c). As can be seen, in this 
frequency range fractons with very different localization lengths are present, in 
qualitative agreement with previous calculations jl^. If, from the modes 
which produce Fig. 3(c), we select those with \B\ < 10~^ and consider the 
distribution of their localization lengths, we obtain Fig. 5(b). Most of these 
low-activity modes are localized on a few masses, and all have practically the 
same frequency (uj « 0.61). Examples of this localized, low Raman active modes 
are in Fig. 6 (a-b) which shows that only few masses are involved in the mode 
(4 and 8 in Fig. 6(a) and 6(b) respectively), while all the other masses are 
practically stationary. In the same frequency range (lu « 0.61) however, there 
are many more modes (see Fig. 5(a)) involving many masses, as exemplified in 
Fig. 6(c). 
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In order to check if a correlation between Raman activity and localization 
length exists, we also studied the D{\B\) relative to groups of modes having 
about the same I. Exception made for modes with / < 3 (which arc those that 
produce the peak in S = of Fig. 3(c)) all the others groups produce gaussian 
distributions for -D(|B|) with exactly the same variance and, therefore, the same 
C{uj). This result was to be expected as a consequence of DILS: in fact, within 
the bond polarizability model, B{p) is the sum of n fluctuating terms {n 
is the number of participating bonds), each having a mean squared amplitude 
« 1/n. This is true if, for a mass involved in mode p, {(ui — Uj)'^) cx (u^), 
because (u^) oc 1/n for normahzed modes. Thus in this case the quantity {B"^) 
does not depend on n or giving rise to the same C{u}) for modes with very 
different localization length. 

Full DID in site percolators was studied as well. What we observe in this case 
is that the longer range of the FD interaction tends to smooth out the extremes 
of almost zero Raman activity observed above, so that it is less simple to single 
out the very localized modes, and the distributions look more like gaussians. 

Bond percolators were also studied. 50 x 50 samples have been used to study 
the low frequency range (Fig. 7(a)). In this case and for FD, as Figs. 7 (b- 
c) show, there are some statistically signiflcant departures (see Table I) from 
the gaussian shape at high frequencies, though the distributions still appear to 
be peaked at zero. This result was expected because in this case there may 
be disconnected pairs of nearest neighbours masses and since connected and 
disconnected nearest neighbours pairs give different DID contributions to \B\, 
a single gaussian distribution is imlikcly to reproduce D{\B\). On the contrary, 
BP produces gaussian-like distributions in this case also. 

4. CONCLUSIONS 

In the present work we studied electrically disordered linear chains with and 
without mechanical disorder, and two dimensional site and bond percolators 
without electrical disorder. 

For the mechanically ordered linear chains, it was analytically shown that 
the distribution _D(|i?|) obtained for each mode by randomly assigning the bare 
polarizability to the masses, is a gaussian centered at zero. Therefore, the 
Raman coupling coefScient of the modes is the variance of the distribution. 
By introducing mechanical disorder, no relevant difference is observed at low 
frequencies, while at higher frequencies i'dSj) shows peaks centered either at 
zero or at finite values of . The latter peaks are due to very localized, Raman 
active modes. Thus, this study singled out two scattering mechanisms: DILS 
(which is intrinsically due to spatial fluctuations) and finite-size scattering. 

The same statistical analysis was performed on site and bond percolators, 
and what emerges is that the DILS mechanism dominates in all frequency ranges 
and for both scattering mechanisms, FD and BP. 

The fiuctuation origin of Raman scattering in these systems is shown by the 
Gaussian shape of the distribution D{\B\). 
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FIGURE CAPTIONS 



Fig. 1. Distribution of the \B\ values relative to 2000 replicas of mechanically 
and electrically disordered chains containing 150 masses each, (a) mi = 1, 
m2 = 2, ai = and a2 = 2 , the first two modes 0.005 < a; < 0.018. (b) mi = 1, 
1712 = 2, ai = 0.9, a2 = 1.1, for all the eigenvalues falling in the frequency range 
0.8994 < w < 0.9050. (c) Same as (b), but for 0.899450 < u)< 0.899475. 

Fig. 2. C(w) for electrically and mechanically disordered linear chains. 
The arrows indicate the frequencies for which the statistical analysis relative to 
figures l(a-b) was performed. 

Fig. 3. Distribution of the l-Baal values computed in the bond polarizability 
model for all modes of 450 replicas of a 20x20 site percolator, (a): 0.08 < a; < 
0.12; (b): 0.37 < w < 0.43; (c): 0.57 < u < 0.63. The arrow in (c) shows the 
number of modes which are left in the first channel neglecting those with \Bxx\ 
and \Byy\ < 10~^. 

Fig. 4. Same as Fig. 3(c), but plotted for \Ba,a\ < 10"^. 

Fig. 5. (a) Distribution of the localization lengths of « 10^ fractons with 
frequency 0.57 < ui < 0.63, i.e. the same which produce Fig. 3(c). (b) Same as 
(a), but considering only the fractons with \Bxx\ and \Byy\ < 10~^. 

Fig. 6. Modes of a 20 x 20 site percolator with (j w 0.61. The localization 
length are: (a) I « 2.1, (b) I « 3.0, (c) I « 7.5. 

Fig. 7. Distribution of \Baa\ values for modes of a bond percolator (FD). 
(a) 50 X 50 percolator in the range 0.014 < to < 0.018. (b) 20 x 20 percolator in 
the range 0.37 < ui < 0.43 (c) 20 x 20 percolator in the range 0.57 < w < 0.63 



TABLE CAPTION 

Table I. Ratios between the first 10 moments, M„, of the distribution D{\B\) 
and the moments of a gaussian having the same zero and second moment of 
D{\B\). The values of the reduced are also reported. 
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